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I.  INTRODUCTION 


The  Computational  Aerodynamics  Branch,  Launch  and  Flight  Division  has  been 
actively  developing  the  capability  to  predict  the  aerodynamics  of  US  Army 
projectiles  using  Computational  Fluid  Dynamics  (CFD)  techniques.  Successful 
prediction  of  Magnus  force  for  a  spinning  artillery  projectile  at  supersonic 
velocities  is  a  prime  example  of  the  powerful  nature  of  these  techniques.1’2 
The  eventual  goal  of  CFD  techniques  Is  to  allow  the  projectile  designer  to 
compute  the  flow  field  about  given  projectile  designs  accurately  and  effi¬ 
ciently  resulting  in  improved  performance  of  future  projectile  designs. 

Currently  under  development  is  the  capability  to  predict  the  supersonic 
aerodynamics  of  finned  missiles  such  as  kinetic  energy  (KE)  penetrators. 3*4 
Kinetic  energy  penetrators  are  long  length-to-diameter  ratio  finned  projec¬ 
tiles  which  are  used  to  defeat  armored  targets  by  punching  holes  through  the 
armor.  Previous  efforts  have  focused  on  computing  the  pitch-plane  aerodynam¬ 
ics  of  these  projectiles.  In  this  report,  several  important  aerodynamic  para¬ 
meters  which  influence  the  roll  characteristics  of  a  fielded  kinetic  energy 
projectile  (M735)  have  been  predicted  using  CFD  techniques.  These  parameters 
include  the  roll  producing  moment  (at  zero  spin  rate),  the  roll  damping 
moment,  and  the  equilibrium  spin  rate,  defined  as  the  spin  rate  for  which  the 
net  roll  moment  is  zero. 

In  contrast  to  conventional  spin-stabilized  shell  configurations,  KE  pene¬ 
trators  rely  on  fins  to  remain  aerodynamical ly  stable,  although  these  projec¬ 
tiles  fly  with  small  spin  rates  to  decrease  the  dispersion  of  the  round.  The 
spin  on  the  projectile  results  from  the  aerodynamic  loads  produced  by  the 
machined  asymmetries  in  the  fin  geometry.  The  geometry  of  the  M735  KE  projec¬ 
tile  is  shown  in  Figures  1  and  2.  The  asymmetries  in  the  fin  geometry  at  the 
leading  and  trailing  edges  of  the  fins  are  typical  of  kinetic  energy  projec¬ 
tiles. 

KE  penetrators  can  be  launched  from  either  smooth-bore  or  rifled  gun 
tubes.  In  the  case  of  the  smooth-bore  gun,  the  projectile  is  launched  with  an 
initial  spin  rate  which  is  nearly  zero  and  rolls  up  to  a  spin  rate  which  is 
nearly  constant  throughout  the  later  stages  of  the  trajectory.  Penetrators 
launched  from  a  rifled  gun  tube  (with  a  slip-ring  sabot)  often  have  an  initial 
spin  rate  higher  than  the  equilibrium  spin  rate  and  roll  down  to  a  nearly 
constant  value.  Aero-ballisticians  describe  the  spin  history  in  terms  of  the 
following  ordinary  differential  equation;5 
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where  p  is  the  spin  rate,  t  is  time,  I  is  the  moment  of  inertia,  and  is  the  i 
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net  aerodynamic  roll  moment  acting  on  the  projectile,  and  pb,  a  ,  M^,  D  and  i 
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diameter,  and  area.  1 

The  net  aerodynamic  roll  moment  is  composed  of  two  components,  the  roll  l 
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asymmetries  in  the  fin  geometry  which  tends  to  induce  spin,  and  the  roll  \ 
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damping  contribution  which  consists  of  pressure  and  viscous  forces  that  oppose 
the  spin.  The  relationship  of  these  contributions  to  the  net  aerodynamic 
roll  moment  is  expressed  below  in  non-dimensional  form. 
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where  C,  is  the  roll  producing  moment  coefficient,  C»  is  the  roll  damping 
o  *p 

moment  coefficient  and  pD/V  is  the  non-dimensional  spin  rate. 

As  is  suggested  by  Equation  (2),  the  roll  producing  moment  can  be  obtained 
by  computing  the  net  roll  moment  at  zero  spin  rate,  while  roll  damping  moment 
is  obtained  by  computing  the  net  roll  moment  on  the  projectile  at  a  fixed  spin 
rate,  and  subtracting  it  from  the  roll  producing  moment  and  dividing  by  the 
spin  rate.  The  equilibrium  spin  rate,  which  occurs  when  the  net  aerodynamic 
roll  moment  is  zero,  is  obtained  by  dividing  the  roll  producing  moment  by  the 
roll  damping  moment. 

In  this  report,  the  net  aerodynamic  roll  moment  is  obtained  by  computing 
the  flow  field  around  the  projectile  using  CFO  techniques.  The  techniques 
applied  in  this  report  are  described  in  the  following  section.  Engineering 
approaches  for  determining  the  roll  characteristics  are  also  briefly  described 
and  are  compared  with  the  computed  results. 

II.  GOVERNING  EQUATIONS  AND  COMPUTATIONAL  TECHNIQUE 

Computation  of  the  viscous  flow  field  about  the  finned  projectile  config¬ 
uration  is  accomplished  by  solving  the  thin-layer  Navier-Stokes  equations 
using  the  Parabolized  Navier-Stokes  technique.  The  computations  have  been 
performed  at  zero  degrees  angle  of  attack  in  a  rotating  coordinate  frame  which 
rotates  at  the  spin  rate  of  the  projectile.  Performing  the  computations  in 
rotating  coordinate  frame  allows  the  time-independent  or  steady-state  govern¬ 
ing  equations  to  be  solved,  at  a  significantly  reduced  computational  cost 
compared  with  computations  performed  in  the  traditional  space  fixed  coordinate 
frame. 

Additional  cost  savings  are  achieved  by  performing  the  computations  in 
cylindrical  coordinates  to  take  advantage  of  the  periodic  symmetry  of  the  cur¬ 
rent  conf iguration  of  interest.  Over  the  cylindrical  portion  of  the  body,  the 
flow  field  is  completely  symmetric  in  the  circumferential  direction,  while 
over  the  finned  portion  of  the  body,  the  solution  is  symmetric  with  a  period 
of  60  degrees,  due  to  the  symmetric  arrangement  of  the  six  fins. 

1.  THIN-LAYER  NAVIER-STOKES  EQUATIONS 


The  set  of  equations  that  describes  three-dimensional  compressible  viscous 
flow  is  referred  to  as  the  Navier-Stokes  equations.  The  Navier-Stokes  equa¬ 
tions  express  the  conservation  of  mass,  momentum,  and  energy  in  the  three 
coordinate  directions.  The  thin-layer  Navier-Stokes  equations,  which  are 
solved  in  this  report,  are  obtained  by  eliminating  from  the  full  Navier-Stokes 
equations,  all  the  viscous  terms  except  for  those  containing  derivatives  in 
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the  direction  nearly  normal  to  the  projectile  body.  For  high  Reynolds  number 
flows  with  no  axial  flow  separation,  the  thin-layer  Navler-Stokes  equations 
are  a  very  good  approximation  to  the  full  Navier-Stokes  equations,  and  can  be 
efficiently  solved  using  available  numerical  algorithms. 

The  Reynolds  averaged  thin-layer  Navier-Stokes  equations,  cast  in 
generalized  cylindrical  coordinates  and  containing  the  rotating  coordinate 
frame  terms,  are  shown  below. 
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Here,  E,  F  and  G  are  the  inviscid  flux  vectors,  Gv  is  the  viscous  flux  vector, 

S  and  Sv  are  the  inviscid  and  viscous  source  terms  due  to  the  cylindrical 

coordinate  formulation,  and  H  is  the  source  term  containing  the  Coriolis  and 
centrifugal  force  terms  which  result  from  the  rotating  coordinate  frame.  Each 
of  these  matrices  are  functions  of  the  dependent  variables  represented  by  the 
vector,  q(p,  pu,  pv,  pw,  e),  where  p  is  the  density,  u,  v,  w,  are  the  velocity 
components  in  the  three  spacial  directions  x,  0,  and  r,  and  e  is  the  total 
energy  per  unit  volume.  Each  of  flux  vectors  and  source  terms  are  defined 
be! ow. 
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where  the  contravariant  velocity  components  are  defined 


U  =  UCX 

V  =  unx  +  vn0/r  +  wnp 
W  =  +  vCg/r  +  wcp 


and  the  metric  quantities  are  defined 
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and  where 
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q2  =  u2  +  v2  +  w2 


Closure  of  this  set  of  equations  is  obtain'd  by  applying  the  ideal  gas 
law,  which  allows  the  pressure  to  be  related  to  the  dependent  variables  in  the 
fol lowing  manner. 


P  =  (y  -  1) [e  -  (p/2)q2] 
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The  turbulent  viscosity,  ut>  which  appears  in  the  viscous  matrices  is 
computed  using  the  Baldwin-Lomax  turbulence  model.6 

2.  THE  PARABOLIZED  NAVIER-STOKES  (PNS)  PROCEDURE 

The  Parabolized  Navier-Stokes  (PNS)  technique  of  Schiff  and  Steger7  is 
widely  used  to  compute  the  supersonic  viscous  flow  about  variety  of  flight 
vehicles  by  solving  the  steady  thin-layer  Navier-Stokes  equations  using  a 
procedure  which  allows  the  solution  to  be  spacially  marched  along  the  body  in 
the  main  flow  direction  due  to  the  parabolic  nature  of  the  governing 


equations.  An  initial  plane  of  data  is  required  to  begin  the  space  marching 
procedure  and  may  be  obtained  either  from  an  auxiliary  calculation  or  from  a 
conical  starting  procedure,  as  has  been  done  for  the  results  presented  here. 
Following  the  approach  of  Schiff  and  Steger,7the  governing  equations,  which 
have  been  modified  to  include  the  Corlolos  and  centrifugal  force  terms,  are 
solved  using  a  conservative,  approximately  factored.  Implicit  finite- 
difference  numerical  algorithm  as  formulated  by  Beam  and  Warming. 8  In  the 
results  presented  here  capturing  of  the  bow  shock  has  been  performed.  Further 
details  of  the  procedure  are  readily  available  In  the  literature. 

For  the  computational  results  presented  here,  the  grid  consisted  of  60 
points  from  the  body  to  the  outer  boundary.  In  the  circumferential  direction, 
the  grlddlng  was  performed  over  a  60  degree  arc  due  to  the  periodic  symmetry 
present  In  the  configuration  of  Interest.  Over  the  axl symmetric  portion  of 
the  body,  six  circumferential  points  were  used.  (This  was  somewhat  arbitrary 
since  the  flow  is  axisymmetric  over  the  body.)  On  the  finned  portion  of  the 
body,  50  circumferential  points  were  used  to  span  the  60  degree  arc. 

Starting  solutions  were  obtained  using  the  conical  step-back  procedure. 
For  each  Mach  number,  starting  solutions  for  the  non-spinning  case  were 
typically  obtained  in  360  seconds  on  a  Cray-2  supercomputer.  Starting  solu¬ 
tions  for  each  additional  spin  rate  were  generated  by  re-converging  the  non- 
spinning  starting  solution  and  typically  required  an  additional  100  seconds. 
Marching  the  solution  over  the  cylindrical  portion  of  the  body  and  the  fins 
required  100  seconds  and  750  seconds,  respectively. 

Solutions  were  obtained  over  the  full  ci rcumferential  plane  (360  degrees) 
to  check  the  implementation  of  the  periodic  boundary  conditions.  The  results 
were  essentially  identical  to  the  periodic  results  and  required  about  five 
times  more  CPU  time.  A  grid  refinement  study  was  also  performed  by  doubling 
the  grid  points  in  the  circumferential  direction.  The  results  showed  only  a 
1.5  percent  variation  in  the  roll  producing  moment  at  Mach  3. 

3.  INVISCID  SPACE-MARCHING  PROCEDURE 

The  viscous  results  have  been  supplemented  by  inviscid  computations.  The 
inviscid  procedure  applied  here  is  based  on  the  PNS  procedure  except  that  the 
viscous  terms  have  been  eliminated,  viscous  sublayer  approximation  ignored, 
and  appropriate  boundary  conditions  have  been  applied  at  the  projectile  sur¬ 
face.  Inviscid  results  have  only  been  obtained  for  zero  spin  rate. 


III.  ENGINEERING  ESTIMATION  OF  THE  ROLL  CHARACTERISTICS 
1.  ESTIMATION  OF  THE  ROLL  PRODUCING  MOMENT 

The  roll  producing  moment  caused  by  the  machined  asymmetries  on  the  fin 
can  he  estimated  using  simple  compressible  flow  theories.  Asymmetries  exist 
at  both  the  leading  and  trailing  edges  of  the  fins,  and  the  contributions  of 
each  to  the  roll  producing  moment  are  discussed  in  the  following  sections. 

a.  Trailing  edge  The  trailing  edge  of  the  fin  is  machined  at  a  constant 
slope  of  seven  degrees  as  shown  in  Figure  2.  This  region  can  be  treated  as  an 
inverted  wedge  and  the  equations  for  Prandtl  Meyer  flow  can  be  applied  to 
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compute  the  pressure  difference  across  the  trailing  edge  bevel,  which  is 
constant  under  the  assumption  of  two-dimensional  flow.  The  roll  producing 
moment  can  be  computed  from  the  pressure  difference  as  shown  below 
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where  A  is  the  area  of  the  trailing  edge  bevel,  Z  is  the  moment  arm  from  the 
projectile  axis  to  the  centroid  of  the  beveled  area,  y  is  the  ratio  of 
specific  heats,  is  the  free  stream  Mach  number,  P,,,  is  the  free  stream 

static  pressure,  AP  is  the  pressure  difference  across  the  beveled  trailing 
edge,  and  N^in$  is  the  number  of  fins. 

By  assuming  free  stream  conditions  in  front  of  the  trailing  edge  bevel 
and  on  the  non-beveled  side  of  the  fin,  the  Prandtl-Meyer  equations  can  be 
applied  to  compute  the  pressure  difference  across  the  fin.  These  relations 
are  available  in  standard  compressible  flow  texts  such  as  Reference  9. 

b.  Leading  edge  The  leading  edge  of  the  fin  for  some  KE  projectiles  may 
have  an  asymmetric  leading  edge  which  can  contribute  to  the  roll  producing 
moment.  Estimates  of  this  contribution  to  the  roll  producing  moment  can  be 
obtained  by  treating  the  leading  edge  chamfer  as  a  two-dimensional  wedge. 
This  component  of  the  roll  producing  moment  can  be  estimated  as  follows, 


lfins»PRZP.  -  UArV,  -  (PL/P.-  DALrL) 


(7) 


where  PR  and  PL  are  the  pressures  on  the  right  and  left  faces  of  the  leading 

edge  chamfer,  AR  and  AL  are  the  areas  of  the  leading  edge  chamfers,  fR  and  VL 

are  the  loc^ions  of  the  centroids  of  the  areas.  The  pressures  PR  and  PL  can 

be  obtained  by  applying  the  oblique  shock  relations,  which  are  also  available 
in  standard  compressible  flow  texts  such  as  Reference  9. 

2.  ESTIMATION  OF  THE  ROLL  DATING  OF  A  KE  FIN 

One  approach  to  estimating  the  roll  damping  of  a  KE  fin  is  to  use  a  strip 
theory  approach,  breaking  the  fin  planform  into  chordwise  strips.  Each  strip 
is  assumed  to  be  a  two-dimensional  flat  plate  at  angle  of  attack,  where  the 
local  angle  of  attack  is  a  function  of  the  local  circumferential  velocity  due 
to  the  spin  and  the  axial  component  of  the  velocity.  The  roll  moment  is  then 
determined  by  integrating  the  lift  on  each  strip  multiplied  by  the  local 
moment  arm.  The  roll  damping  moment,  which  is  the  variation  in  the  roll 
moment  with  spin,  can  be  determined  by  dividing  the  roll  moment  by  the  non- 
dimensional  spin  rate,  since  roll  moment  on  the  flat  plate  fin  will  be  zero  at 
zero  spin  rate.  This  assumes  that  the  beveled  trailing  edge  has  little  effect 
on  the  roll  damping.  This  is  consistent  with  the  linearized  potential  theory 
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applied  In  the  context  of  the  spinning  KE  penetrator,  which  separates  the 
trailing  edge  deflection  effect  (roll  producing)  from  the  angle  of  attack 
effects  (roll  damping). 

From  linearized  potential  theory,9  the  lift  coefficient  on  a  two- 
dimensional  flat  plate  at  angle  of  attack  may  be  expressed  as  follows 


C 


L 


(8) 


By  assuming  that  the  flow  conditions  in  front  of  the  fins  are  approximate¬ 
ly  free  stream,  the  lift  on  the  strip  can  be  written. 


L  =  _Mr_L  I  pjfcajfrr) 

/M2  .  1  2 

op 


(9) 


where  c(r)  and  a(r)  are  the  local  chord  length  and  angle  of  attack  respec¬ 
tively. 

The  roll  moment  coefficient  is  obtained  by  integrating  the  product  of  the 
lift  and  the  moment  arm,  and  non-dimensional izing  appropriately. 


The  local  angle  of  attack  can  be  approximated  as  follows. 


(11) 


The  roll  damping  coefficient  is  obtained  by  dividing  the  roll  moment 
coefficient  by  the  non-dimensional  spin  rate. 


IV.  RESULTS 

Computations  have  been  performed  to  determine  the  following  aerodynamic 
parameters  which  influence  the  roll  characteristics  of  kinetic  energy  penetra- 
tors:  the  roll  producing  moment;  the  roll  damping  moment;  and  the  equilibrium 
spin  rate.  Viscous  computations  were  performed  over  a  range  of  Mach  numbers 
(M  =  3.0  to  4.5)  and  spin  rates  (pO/V  =  0.  to  pO/V  =  0.0964)  for  free-flight 
(sea-level)  atmospheric  conditions.  Inviscid  computations  were  performed  to 
supplement  the  viscous  computational  results  at  zero  spin  rate.  Unless  other¬ 
wise  noted,  the  computational  results  discussed  below  refer  to  those  obtained 
using  the  viscous  PNS  procedure. 
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1.  ROLL  PRODUCING  MOMENT 


Figure  3  shows  the  development  of  the  computed  roll  producing  moment  coef¬ 
ficient  as  a  function  of  the  axial  position  along  the  body  at  Mach  3  and  zero 
spin  rate.  The  computational  result  shows  that  for  the  W735  geometry,  the 
roll  producing  moment  is  produced  almost  entirely  by  the  fin  asymmetries  at 
the  fin  trailing  edge  and  that  the  fin  asymmetries  at  the  fin  leading  edge 
contribute  very  little  towards  producing  a  roll  moment.  This  is  consistent 
with  the  estimates  of  the  roll  producing  moment  obtained  using  the  engineering 
approach. 

Figure  4  shows  the  roll  producing  moment  coefficient  as  a  function  of  Mach 
number  as  predicted  by  both  viscous  and  in viscid  computation  and  by  the  engi¬ 
neering  approach.  The  trend  predicted  by  the  computational  and  engineering 
approaches  compare  well,  though  the  engineering  approach  predicts  a  roll 
moment  30  to  50  percent  greater  than  in  the  viscous  computation.  The  roll 
moment  predicted  by  the  inviscid  computation  ranges  from  15  to  28  percent 
greater  than  the  viscous  result. 

Figures  5a  and  5b  display  the  pressure  difference  across  the  fin  (from  one 
side  to  the  other)  in  the  vicinity  of  the  trailing  edge  as  a  function  of  the 
distance  from  the  leading  edge.  This  pressure  difference  is  the  mechanism 
that  causes  the  roll  producing  moment.  Pressure  profiles  at  seven  span  loca¬ 
tions  are  shown.  Section  1  is  located  at  one-eighth  of  the  distance  between 
the  fin  root  and  the  fin  tip,  while  section  7  is  located  at  seven-eighths  of 
this  distance.  The  other  five  sections  are  spaced  equally  in-between  as  shown 
in  Figure  6.  Figures  5a  and  5b  display  the  inviscid  and  viscous  results, 
respectively.  A  reference  value  corresponding  to  the  pressure  drop  across  a 
two-dimensional  wedge  is  also  plotted  on  both  of  these  figures. 

Both  sets  of  results  show  that  in  front  of  the  trailing  edge  chamfer,  the 
magnitude  of  the  pressure  difference  is  nearly  zero,  but  rises  sharply  on  the 
trailing  edge  chamfer.  The  inviscid  results  attain  the  reference  value  except 
at  positions  near  the  tip  of  the  fin.  The  viscous  results  approach  the  refer¬ 
ence  value  in  a  more  asymptotic  manner  due  to  the  presence  of  the  boundary 
layer.  In  the  viscous  results,  tip  effects  also  produce  a  smaller  pressure 
differential  at  section  7.  The  viscous  results  also  show  a  slightly  smaller 
pressure  differential  at  section  1.  This  can  be  attributed  to  the  thick 
boundary  layer  on  the  cylindrical  portion  of  the  body.  The  thickness  of  this 
boundary  layer  can  be  seen  in  Figure  7,  which  shows  the  Mach  number  contours 
at  an  axial  station  in  the  vicinity  of  the  trailing  edge. 

These  detailed  pressure  profiles  help  explain  the  differences  between  roll 
producing  moments  predicted  by  the  viscous,  inviscid,  and  engineering  esti¬ 
mates.  The  inviscid  roll  producing  moment  is  larger  than  the  viscous  roll 
producing  moment  because  the  effect  of  the  body  and  fin  boundary  layers  on  the 
driving  pressure  difference  is  neglected,  due  to  the  inviscid  assumption.  The 
engineering  estimate,  which  assumes  that  the  pressure  difference  across  the 
trailing  edge  chamfer  is  the  reference  value,  over-estimates  the  roll  produc¬ 
ing  moment  because  tip  effects  as  well  as  viscous  effects  are  completely 
ignored. 


'.V 


si 


2.  ROLL  MOMENT  ON  A  SPINNING  PROJECTILE 


For  each  Mach  number,  computation  of  the  flow  field  over  the  bod[y  has  been 
performed  for  a  number  of  fixed  spin  rates.  The  direction  of  spin  is  in  the 
same  direction  as  the  zero- spin  roll  moment. 

Figure  8  displays  the  development  of  the  net  roll  moment  coefficient  over 
the  projectile  at  a  flight  velocity  of  Mach  3  and  rolling  at  a  non-dimensional 
spin  rate  of  .008568  (40  rps).  This  figure  shows  that  relatively  little  roll 
damping  is  produced  by  the  axisymmetrlc  portion  of  the  body.  Once  on  the  fin, 
the  net  roll  moment  shows  an  Initial  rapid  decrease  as  a  result  of  the  roll 
damping  contribution  which  opposes  the  spin.  Near  the  aft  end  of  the  projec¬ 
tile,  the  net  roll  moment  begins  to  Increase  as  the  roll  producing  moment 
predominates  over  the  roll  damping  contribution.  For  the  particular  spin  rate 
shown  here,  the  total  net  roll  moment  Is  positive.  Indicating  that  the  equili¬ 
brium  spin  rate  will  be  greater  than  this  particular  spin  rate. 

The  computed  net  roll  moment  coefficient  as  a  function  of  spin  rate  at 
Mach  3.0  is  shown  In  Figure  9.  Clearly  shown  is  the  decrease  In  net  roll 
moment  with  spin  rate.  This  decrease  is  due  to  the  increase  in  the  roll 
damping  contribution  of  the  net  roll  moment  as  the  spin  rate  increases.  The 
results  confirm  that  net  roll  moment  varies  linearly  with  spin  rate,  as 
discussed  previously. 

The  slope  of  curve  shown  in  Figure  9  represents  the  roll  damping  for  this 
Mach  number.  Similar  computations  have  been  performed  at  the  other  Mach 
numbers  and  the  roll  damping  computed  for  each  of  these  Mach  numbers.  The 
computed  roll  damping  coefficient  as  a  function  of  Mach  number  is  shown  in 
Figure  10.  Also  shown  Is  the  roll  damping  coefficient  as  predicted  using  the 
engineering  estimate.  The  engineering  estimate  predicts  a  roll  damping  coef¬ 
ficient  which  is  28  to  34  percent  greater  than  the  computation,  though  the 
trend  of  both  predictions  is  consistent. 

Figure  9  also  shows  that  at  Mach  3,  the  equilibrium  spin  rate,  (the  spin 
rate  at  which  the  net  roll  moment  becomes  zero)  occurs  at  a  non-dimensional 
spin  rate  of  «onor>8,  corresponding  to  42.4  revolutions  per  second.  The  non- 
dimensional  equilibrium  spin  rate  as  a  function  of  Mach  number  is  shown  in 
Figure  11.  Equilibrium  spin  rates  obtained  from  the  engineering  approach  are 
also  shown  in  this  figure.  Differences  of  up  to  seventeen  percent  between  the 
two  approaches  are  observed.  The  computational  results  show  less  variation  of 
the  dimensional  equilibrium  spin  rate  with  Mach  number  than  does  the  engineer¬ 
ing  estimation,  as  shown  in  Figure  12. 

3.  TRAJECTORY  SIMULATION 

The  roll  history  of  the  projectile  has  been  computed  by  solving  Equations 
(1)  and  (2).  The  effect  of  the  projectile  deceleration  due  to  drag  (65  km/sec 
per  kilometer)  has  been  considered.  Roll  histories  have  been  obtained  using 
both  the  computed  roll  coefficients  and  the  coefficients  obtained  using  the 
engineering  approach.  The  roll  histories  were  obtained  for  a  launch  velocity 
of  Mach  4.25,  atmospheric  conditions,  and  zero  Initial  roll  rate.  The 
results,  shown  in  Figure  13,  show  a  rapid  spin-up  over  the  first  half  kilo¬ 
meter  of  flight,  followed  by  a  slow  decrease  in  spin  rate  over  the  rest  of  the 
trajectory.  Both  trajectories  are  similar,  though  the  results  obtained  using 
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the  engineering  estimated  coefficients  show  a  slightly  higher  spin  rate 
throughout  the  trajectory.  The  results  also  show  that,  after  about  0.75  kilo¬ 
meters,  the  projectile  is  flying  at  the  equilibrium  spin  rate  where  the  spin 
rate  is  just  a  function  of  the  local  Mach  number.  This  is  shown  in  Figure  14, 
where  the  equilibrium  spin  rate  and  the  spin  rate  from  the  trajectory  simula¬ 
tion  as  a  function  of  Mach  number  are  displayed.  The  results  show  that  after 
the  initial  spin-up  transient,  the  projectile  is  flying  at  the  equilibrium 
spin  rate. 

V.  CONCLUSION 

The  roll  characteristics  (roll  producing  moment,  roll  damping  moment,  and 
equilibrium  spin  rate)  of  a  kinetic  energy  penetrator  have  been  obtained  using 
Computational  Fluid  Dynamics  techniques.  The  computed  results  show  that  vis¬ 
cous  computation  is  required  to  compute  the  roll  characteristics  accurately, 
based  on  the  differences  between  the  viscous  and  Inviscid  results. 

The  computational  results  have  been  used  to  benchmark  engineering 
approaches  for  estimating  these  roll  characteristics.  For  the  particular 
projectile  examined  here,  the  engineering  approaches  over-estimated  the  roll 
producing  moment  by  up  to  50  percent  and  the  roll  damping  moment  by  as  much  as 
34  percent,  though  the  trend  predicted  by  computation  and  engineering  estima¬ 
tion  is  similar.  The  predicted  values  of  equilibrium  spin  rate  differ  by  a 
maximum  of  17  percent,  though  the  computed  results  show  a  greater  variation 
with  Mach  number  than  do  the  values  determined  from  the  engineering  estimation 
approach. 
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